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Abstract. An overview is given about recent developments and results of comprehensive simu- 
lations of magneto-convective processes in the near-surface layers and photosphere of the Sun. 
Simulations now cover a wide range of phenomena, from whole active regions, over individual 
sunspots and pores, magnetic flux concentrations and vortices in intergranular lanes, down to the 
intricate mixed-polarity structure of the magnetic field generated by small-scale dynamo action. 
The simulations in concert with high-resolution observations have provided breakthroughs in our 
understanding of the structure and dynamics of the magnetic fields in the solar photosphere. 
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1. Introduction 

Solar magneto-convection is a rather extended area of research. It covers topics ranging 
from the generation of the large-scale magnetic field by the solar dynamo in the deep 
convection zone to the excitation, propagation, and dissipation of currents and MHD 
waves in the chromosphere and beyond. Obviously, this wide field cannot be reasonably 
covered in a sin gle review. R ecent reviews of various asp e cts of solar magne t o-con vection 
can be found in(Fan| (|2009l ) . ICharbonneaul (fioToh . ISteinl (|2012h . and lWeisd (|2012h . 

Here I shall restrict the discussion to the solar photosphere, the region where the 
best understanding of the magneto-convective processes has been achieved so far. The 
significant progress was made possible through the interplay between high-resolution 
observations, idealized theoretical models, and comprehensive numerical simulations. The 
latter provide a 3D view of the processes underlying the observations. Simulations can be 
validated (albeit not in a strict sense) by detailed comparison of synthetic observational 
quantities based on simulation results with real observations. At the same time, this 
comparison provides a firm basis for the interpretation of the observations in terms 
physical quantities and processes. 

I prefer to qualify these simulations as 'comprehensive' in favor of the often-used term 
'realistic', the reason being that the values of important non-dimensional numbers (such 
as the kinetic and magnetic Reynolds numbers, or the magnetic Prandtl number) that 
can be reached in the simulations differ from the realistic values by (many) orders of 
magnitude. Typically, some kind of hypcrdiffusion or subgrid modelling is used to ac- 
count for the unresolved spatial scales, but in the absence of a theory of turbulence the 
validity of such quasi-heuristic approaches cannot be strictly proven. On the other hand, 
the term 'comprehensive' expresses the fact that these simulations aim at representing all 
relevant physical processes and conditions in the photosphere, i.e., compressible MHD in 
3D domain, partial ionization and molecule formation, as well as non-grey and non-local 
radiative transfer. A number of numerical c odes are used to obta in comprehensive sim- 
ulations, amone 



(iFrevtag et al.l 



g the m the ANTARES co de (IMuthsam et alJl2010h. t he C05BOLD code 
l2012h. the MURaM code (<Voglerll2003 HVogler et af ] l2005l). the PENCIL 



code (jBrandenburg fc Doblerl 120021). the STA GGER/BIF ROST code (GudikscnetaJ 



2011 ). as well as the codes of Robinson et al.l (|2003l ). and I Jacoutot et aL ( 2008 ). There 
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are also codes with rather coarse tre atment of radiation, such as the RADMHD code 
(lAbbett 2007 ; Abbett fc Fisher 2012 ) with a local radiative cooling term or the code of 
Ustvugovl (|2010l ). which uses the diffusion approximation. Unfortunately, not all groups 



have given detailed descriptions of the numerical methods and approximations used in 
their codes and shown test results in the literature, so that it is difficult to judge and trust 
the validity of their simulation results. On the other hand, the STAGGER, C05BOLD 
and MURaM groups groups have embar ked on a cross-va lidation of their codes by a 



detailed comparison of numerical results (jBeeck et al.ll2012l ) 



Even when restricting oneself to the photosphere, a complete review of observational, 
theoretical, and simulation results is far beyond the scope of this contribution. Limitations 
of space (and human time) force me to restrict myself mainly to the discussion of recent 
simulation results, with very limited references to observation and theoretical concepts. 
Furthermore, I am afraid that this text is somewhat biased towards results obtained with 
the MURaM code - not because this code would be better than others, but just because I 
am most familiar with its results. What follows is organized according to the spatial scale 
of the phenomenon: from whole active regions down to the fine threads of the magnetic 
field generated by small-scale dynamo action. 



2. Active regions 

In the past decade, the computational resources available for comprehensive simu- 
lations have increased dramatically. It has become possible to simulate the emergence 
and development of whole acti ve regions with suns pots, pores, and plage areas. This 
line of research was started by Ch eung et al.1 (|2007l ). who considered the emergence of 
a twisted flux tube carrying a longitudinal magnetic flux of 10 19 Mx. The initially co- 
herent flux tube becomes undulated and fragmented by the vigorous convective flows, 
so that the flux emergence takes place in the form of numerous patches of horizontal 
field connecting small bipoles in a 'salt-and-pepper' pattern. The disturbed granulation 
at the emergence site shows elongated gran ules threaded by dark lanes. Many of these 
featu res are also found in observati ons (e.g.. ICheung et al. 20081 Yelles Chao u che et al. 
200 j. Similar results were found by IStein et al.l (|201ll ). iFang et al.l (|2010l l2012h . andbv 



Martinez- Svkora et al. ( 2008 . 2009t ). The latter authors also considered the evolution of 



t he coronal field durin g flux emergence and the development of spicule-like structures. 

Cheung et al. I (l2010h simulated the emergence of a flux tube with 7.6 x 10 21 Mx. As the 
emergence proceeds, the pattern of small bipoles at the visible solar surface increasingly 
reorganizes until the magnetic flux distribution reflects the large-scale bipolar structure 
of the underlying Omega-shaped flux tube. At the same time, a pair of sunspots forms 
through the coalescence of pores and smaller flux fragments, so that the appearance 
of a full-fledged active region develops. These simulations also solved the long-standing 
problem of how the enormous amount of mass contained in a deep-lying flux tube is 
removed as it reaches the low-density regions near the surface. It turns out that much 
of this mass is accumulated in sinking U-shaped loops that reconnect near the surface 
and thus release the mass from the field line s of the surface field. Observational evidence 
for this process has recently been reported (|Centenoll2012l) . Simulations by M. Cheung 
and M. Rempel (not yet published) show that an initial twist is not essential for the 
formation of coherent sunspots from a flux tube inserted at a depth of about 16 Mm. 

While the simulations discussed so far assumed a large-scale flu x tube to coherently 
enter the simulation domain from below, IStein fc Nordlundl (l2012h assumed horizontal 
magnetic field to be advected by upflows at the bottom boundary of the computational 
domain. They found the development of a bipolar structure at the surface that resembled 
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a small active region with a few pores. Although this simulation shows that a bipolar 
region can form without inserting a coherent flux tube beforehand, some caution is indi- 
cated when interpreting this result: the horizontal field advected into the computational 
box has a fixed direction and strength, so that the simulation implicitely assumes the ex- 
istence of an infinite reservoir of strictly organized and coherent horizontal magnetic flux 
below the computational domain, in fact something like a huge homogeneous flux sheet. 
Consequently, the orientation of the developing bipolar region is along the direction of 
the assumed horizontal field below the simulation box. 



3. Sunspots and pores 

Comprehensive simulations of sunspots have seen tremendous progress in the last few 
years. Already studies of 'sunspot slabs' at ra ther coarse resolution and limited box size 
(jHeinemann et al.ll2007t iRempel et al.ll2009bl ) provided crucial insights in the processes 
responsible for the pen umbral filamentation and the Evershed effect. The first simula- 
tions of full sunspots by Rempel et al.l (|2009al ) and their further development and detailed 
analysis bv lRempell (|2011aUbl 12012^ revealed the nature of the magneto-convective pro- 
cesses leading to the characteristic structure of sunspots. While the bright umbral dots 
represent hot upfl ow plumes, whose lateral expansion is strong ly constricted by the ver- 



tical umbral field (Schiisslcr & Voglcr 2006; Bharti et al. 20f0), the inclined field in the 



penumbra provides a path for the unimpeded expansion of hot upflowing material: the 
direction along the field. As a consequence, sheet-like convective upflows in the penumbra 
appear as bright, elongated filaments with strong outflows (representing the Evershed ef- 
fect), which feed downflow patches in the outer part of the penumbra. Part of the rising 
plasma turns over laterally (perpendicular to the filament axis) and descends in downflow 
lanes next to the bright filaments. The expanding outflow stretches and expels the mag- 
netic field, so that the outflow channels exhibit relatively weak (albeit non-vanishing) and 
almost horiz ontal magnetic field. All these properties are consistent with observational 
results (e.g..lBqrrero fc IchimotolboilL IRempel fc Schl chenmaierll201ll : ljoshi et al.ll20lH 



S charmer et alj 



2011). 



It is observationally well established that the formation of big sunspots invariably is 
connected to major events of magnetic flux emergence. This is also reflected in simula- 
tions of sunspot formation and structure, which always require either the advection of 
well-organized magnetic flux into the computational domain (as in the simulations of 
active-region formation discussed in Sec. [2]) or the assumption of a sufficiently big con- 
centration of vertical magnetic flux as initial condition for the simulation. The situation 
is different with pores, small dark flux concentrations without penumbrae and mostly 
not exceeding the area of a few granules. Such structures are observed to form and decay 
quasi-spontaneously in mature plage areas, without any associated emergence of new 
magnetic flux. Magneto-convection simulations in sufficiently deep boxes also show the 
formation of pore-like dark structures if enough magnetic flux is available, either in the 



form of a pres et vertical flux through the computational box ( Kitiashvili et al. 20101 : 



Ustvugovll2010l) or by the perpetual advection of horizontal magnetic flux through the 
bottom boundary ([Stein et al.l 120111 ). Previous simulations in shallow boxes extending 
to a few Mm depth below the visible solar surface did not show the formation of dark 



structures, apart from occasio nal small 'micropores' (e.g., Bercik et al.ll2003 : Vogler et al 



20051: IStein fc Nordlundl[2006h . The reason why shallow simulations miss the formation of 



proper pores is that it requires horizontal flows on a sufficiently big spatial scale in order 
to collect the necessary sufficient amount of magnetic flux. Since the horizontal scale of 
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Figure 1. Maps of brightness (bolometric intensity, left panel) and vertical magnetic field at 
the optical solar surface (right panel) from a snapshot of a magneto-convection simulation in a 
49 x 49 Mm 2 wide and 15 Mm deep box permeated by a horizontally averaged vertical field of 
100 G. 



the convective flows increases with depth (see, e.g., Fig. 6 in lNordlund et al.ll2009l ) . only 
sufhcently deep boxes (at least 6 Mm) provide the conditions for pore formation. 

A word of caution is in order here. The maximum horizontal scale of the convective 
motions seen in simulations is roughly equal to the depth. This means that the horizon- 
tal extension of the computational box should be at least 3-4 times its depth in order 
to provide enough volume for a development of a large-scale flow pattern that is not 
restrained by the periodic bounday conditions at the side walls of the box. Simulation in 
boxes with smaller aspect ratios most probably provide unrealistic results. 

What determines the size (or amount of flux) of pores and why do big sunspots not 
form in this way? This depends basically on the amount of magnetic flux available (mean 
flux density) and on the relation of the advection time of the horizontal flow (horizontal 
scale divided by flow speed) and the lifetime of the flux concentration: the advected 
flux per unit time has to be at least equal to the flux lost by 'turbulent erosion' of 
the flux concentration ( Cameron et al.l [20071) . As a result, for the mean vertical flux 
densities found in plages (a few hundred G), the biggest flux concentrations formed by 
flux advection do not exceed the size of a few granules. As an example, Fig.[T]shows maps 
of the emergent intensity (left panel) and of the surface distribution of vertical magnetic 
field (right panel) from a simulation (run by M.C.M. Cheung) in a box of 49 x 49 Mm 2 
horizontal extension reaching down to 15 Mm depth, so that nearly supergranular-size 
flows can develop. The average vertical field strength is 100 G, thus representing a typical 
plage region. The simulation formed a number of pores, which qualitatively resemble those 
occuring in plage areas. Furthermore, the distribution of the magnetic field shows a multi- 
cellular pattern, reflecting the various horizontal scales of the convective flow. The biggest 
cells have a scale of about 15 Mm, corresponding to the depth of the computational box. 
This pattern is remarkably similar to those seen in plage observations. This result opens 
the possibility to obtain information about the flows in the upper part of the convection 
zone by surface observations and comparison with simulations, providing a supplement 
to hclioseismic inversions. 

Fig. [2] shows the distribution of the vertical magnetic field component on a vertical 
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Figure 2. Map of the vertical magnetic field component, B z , in a vertical cut along a horizontal 
line at x — 27.5 Mm in Fig. Q] (right panel). Reddish colors indicate positive values of B z , bluish 
colors negative values. The zero level of the depth scale corresponds to the optical surface (mean 
level of optical depth unity). 

section of the simulation box. The section cuts through surface flux concentrations of var- 
ious sizes. While the small-scale flux concentrations in intergranular lanes are restricted 
to the near-surface layers, bigger flux concentrations extend deeper, reflecting the depth 
of the flow pattern that accumulated them in the first place. The big magnetic flux 
concentration corresponding to the pore near the right-hand edge of the plot stretches 
down almost to the bottom of the computational domain. While the magnetic field is 
rather organized in the surface layers, the deeper parts of the box show patches with 
a chaotic pattern of small-scale field structures, resulting from the turbulent nature of 
the convective downflow regions. These are the sites of small-scale dynamo action (see 
Sec,©. 



4. Small-scale flux concentrations 

Starting with the pioneering work of lNordlundl (Il983l Il985bl) . comprehensive 3D sim- 
ulations have revealed the processes responsible for the concentration of most of the 

'non-turbulent' part of the magnetic flux at the solar surface into small-scale structures 

I II 1 

with kG field strength located in intergranular downflow lanes (e.g., iSteiner et all 1996. 
19981: iBercik et al.lll998l: IStein et alj|2002t IVogler fc Schussler ll200a IVogler et al-l booa 
Schaff enberger et al. 2006 ; Stein fe Nordlundl 20061 ) . The simulations confirmed the valid- 



ity of the basic theoretical explanations for the intermittent structure of the solar surface 
field, which were developed already in the 60s and 70s of the last century. The first con- 
ce pt is flux expulsion by c onvective flow s in an electrically conducting fluid, pioneered 
bv lWeissl (Il964l Il966l) and lParker! (|l963lh The advection of magnetic flux by horizontal 
flows and its concentration in convective downflow areas concentrates the field roughly up 
to equipartition field strength, for which the magnetic energy density equals the kinetic 
energy density of the flow. The concept of convective collapse is based in the strongly 
superadiabatic startification of the subsurface layers, which renders equipartition flux 
concentrations unstable with respect to downflow along the field lines ( Webb fc Roberts! 
1978USpruit fc Zweibelllil)79ri . This evacuates the upper part of the flux concentrations, 



so that the gas pressure of the s urrounding p l asma compr esses the flux concentration 
until kG field strength is reached ([Parker Il978l ISpruitl[l~979( ) . The depression of the sur- 
face of optical depth unity in the evacuated flux concentration leads to la teral heating 
of its interior by th e hot walls and its appearance as a bright structure (jSpruitl 11976 ; 
Deinzer et al.l Il984l) . Quantitative comparison of observations and comprehensive nu- 
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mcrical simul ations confirmed the relation between kG magnetic flux concentrations and 



bright points (Schussler ct al. 2003; Shelyag ct al. 2004: Lagg ct al. 


2010; 


Rohrbein et al. 


2011; Danilovic et al. 


2012), faculae (Keller et al. 2004; Carlsson et al. 


2004), and the 


SDectro-Dolarimetric signatures of these structures (Shelvag et al.ll2007l). 



5. Vortices 
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Figure 3. Volume rendering of the swirling strength (for a definition, see Moll et al. 2011al ) for 
a non-magnetic simulation (upper panel) and for a simulation with an average vertical magnetic 
field of 200 G (lower panel), viewed from one side of the computational box. Shown is the upper 
half (6 x 6 x 0.8 Mm 3 ) of the box, the bottom plane corresponding roughly to the optical surface. 
While the vertical vortex tubes turn over near the surface in the non-magnetic case, they align 
with the vertical flux concentration and protrude high into the atmosphere when a background 
magnetic field is present. 



The formal equivalence of the vorticity equation and the induction equation suggests 
that vertical vorticity should be expelled and concentrated in intergranular downflows 
in similar ways as magnetic flux. It is no surprise, therefore, that vort ical features and 
whirl flows are regularly found in comprehensive convection simulati ons (jNor dlundll985at 
Stein fc Nordlundlll998t iMoll et ahll2011al IKitiashvili et "al"]|2012al ). While most of the 



strong vertical vortices have diameters below 100 km, probably too small to be detected 
with existing telecopes, larger-scale v o rtex f l ows have been reported by observers (e.g. , 
" Bonet et al] 120081 l20ld lAttie et al.l 120091; IVargas Dommguez et al" 



Brandt et al. 1988 



20111) . In addition. ISteiner et al.l (|2010l ) detected horizontal vortices in granules, both in 



simulations and in observations with the Sunrise balloon telescope. 

If a sufficient amount of vertical background flux is prese nt, the vortices are closely 
associated with magnetic flux concentrations. IVogler (2004) first found 'magnetic vor- 
tices' in magneto-convection simulations. Recently, this topic found consid erable interest, 
also in connection with chromospheric and coronal heating p rocesse s (e.g., Shelvag et al 



201 "it IKitiashvili et al.ll2012bHSteiner fc RezaeHl2012HWedemever-Bohm et al.ll2012l) . The 



simulations of lMoll et ahl (|2012l)"" indicate that the dynamical structure of the upper pho- 
tosphere and chromosphere of 'quiet' regions (i.e., areas with a mean signed field below 
about 20 G) is significantly different from that of more strongly magnetized regions: 
while th e former are dominated by patterns of moving shock fronts (as noted previ- 
ously by IWedemever et al. 1 12004 IWedemever-Bohm1l2010l ). the more magnetized areas 
show numerous vertical vortices, which are strongly associated with the small-scale mag- 
netic flux concentrations in intergranular lanes. As illustrated in Fig. [31 vertically orien- 
tated vortices in the surface layers of weakly magnetized regions do not extend high into 
the atmosphere, b ut turn over near the optical surface and form low-lying vortex loops 



(jMoll et al.ll2011al ). In more strongly magnetized regions (network, plage), the vertical 
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vortices align with the magnetic flux concentrations and reach high into the atmosphere, 
where the dissipation of kinetic energy leads to intense local heating. Fig. [4] illustrates the 
different dynamical character of the upper photosphere/lower chromosphere in weakly 
and strongly magnetized regions: while shocks dominate the former, hot vortices provide 
the governing structure in the latter case. 




Figure 4. Horizontal cuts (size 6 Mm x 6 Mm) through the upper photospheric layers of a 
nonmagnetic simulation (upper row) and a simulation with an average vertical magnetic field of 
200 G (lower row). Shown are maps of temperature (left panels), ranging from about 3000 K to 
over 7000 K, and of horizontal flow speed (middle panels), which reaches values up to f 5 km-s -1 ). 
The upper right panel shows a (bolometric) brightness map for the hydrodynamic case and the 
lower right panel gives a map of the vertical magnetic field in the upper photosphere of the 
magnetic run. The hydrodynamic case is dominated by shock fronts while vortices prevail in the 
magnetic simulation. Both kinds of structures show strong local heating of the plasma. 



6. Small-scale dynamo 

There are a number of observational indications for the existence of a considerable 
amount of mag netic flux in the form of a small-scale turbu lent field in the solar near- 



surface layers ([Sanchez Almeida fc Martinez Gonzalez! 120111 see also the contribution 



of S. Tsuneta in this volume). This compon ent of the solar surface flu x seesm to be 
independent of the 11-year solar cycle (e.g., iTruiillo Bueno et al. 2004), indicating a 
different source mechanism. Recently, D. Biihler (in preparation) systematically analyzed 
Hinode/SP data for the period between 2007 and 2012, which covers the extended solar 
minimum and the rise of the current cycle. He found that the area fraction of weak, but 
significant, linear and circular polarization signals in quiet internetwork regions did not 
change during this time. This suggests that the turbulent magnetic field is not connected 
to the solar cycle, but possibly generated by small-scale dynamo (SSD) action. 

Characteristic for a SSD driven by a turbulent flow in an electrically conducting fluid 
is that the magnetic field is generated at spatial scales that are much smaller than that of 
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the energy-carrying eddies of the flow (i.e., the integral scale of the turbulence). Although 
there is some dependence on the value of the magnetic Prandtl number (ratio of kinematic 
viscosity to magnetic diffusivity) , SSD action p robably occurs in all tur bulent flows of 
sufficiently high magnetic Reynolds number (see Brandenburg et al]|2012l and references 
therein) . While direct numerical simulations demonstrated SSD action in various set tings 
since the 1980s (e.g., iMeneguzzi et al.l Il98ll: ICattaneol fl999t iBushbv et all l2012h . the 
effect was also found in large-eddy simulations (Bo ussinesq or anelastic) in spherical 



shells carried out to mo del the solar convection zone (jGilman fe Miller 1981 : Glatzmaier 
198.4 iBrun et~aT1l2004h . 

Comprehensive simulations of solar near-surface convection indicate that the observed 
turbulent field could indeed be a product of a SSD action driven by the turbulent inter- 
granul ar downflows ( Vogler fc Schussler 2007 : Pietarila Graham et al. 2010t Moll et al 



2011blh The characteristic properties of a magnetic field generated by a SSD can explain 



the o bserved strong horizontal fields observed in the middle photo sphere ([S chussle r fc Vogler 



2008) and the weak signals d etected with sensitive polarimeters ([Pietarila Graham et al 



2009t iDanilovic et al.l[2010al lbl). While most of the field due to the SSD is of mixed polar 



3gl' 
J. 



ity on small scales and has a strength of some tens to a few hundred Gauss, occasionally 
enough unipolar flux is being assembled by the granular flows to form a kG flux concen- 
tration appearing as a bright point in the visible light. Therefore, observations of such 
features in quiet internetwork areas on the Sun are consistent with SSD action. 




Figure 5. Snapshot from a simulation of small-scale dynamo action in the solar near-surface 
layers. The computational box is 12 x 12 Mm 2 wide and 6.1 Mm deep. The panels show horizon- 
tal cuts of the vertical field component (upper panels; black and white indicates negative and 
positive polarity, respectively) and of the vertical flow velocity (lower panels; light shades indi- 
cates upflows, dark shades downflows). The cuts were taken at a depth of ~ 5 Mm (left panels) 
and 2.5 Mm (middle panels) below the average height of the optical surface, as well as at the 
optical surface (right panels). The dynamo-generated field is associated with the downflows in 
the deeper parts of the domain and thus exhibits a 'mesogranular' pattern at the surface. 



Fig. [5] shows a snapshot from a dynamo simulation carried out in a deeper and wider 
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computational box than previous simulations. The maps show a close association between 
the dynamo-generated field, which exhibits the characteristic mixed-polarity pattern of 
small-scale dynamo action, and downflow areas. This is due to the fact that the small-scale 
dynamo mainly works in the turbulent downflows. The long-lived, large-scale horizontal 
convective flow pattern in the deep layers is r eflected in the 'mesogranula r' distribution 
of the magnetic field in the surface layers (cf. lYelles Chaouche et al. 2011 ) . The scale of 
this pattern is determined by the relation between the horiontal advection time and the 
timescale of flux generation (growth rate of the SSD) in the intergranular downflows. 



7. Concluding remarks and outlook 

The combined efforts of high-resolution observations and comprehensive numerical 
simulations have tremendously improved our understanding of the magneto-convective 
processes in the near-surface layers and the photosphere of the Sun. The identification of 
key physical processes through the simulations and the detailed quantitative comparison 
between observational data and numerical results was the key factor for this progress. It 
should be noted that such a comparison needs to be done in a proper way, which means 
that observational data (e.g., intensity images or Stokes profiles) should be compared with 
the corresponding quantities derived from the simulation results (i.e., synthetic Stokes 
profiles) after taking into account instrumental effects (e.g., convolving with point spread 
functions, adding straylight and noise, rebinning to detector pixel size). Comparing sec- 
ondary observational products such as inversion results directly with physical quantities 
from a simulation can lead to misleading conclusions. 

Where do we go from here, which are the frontiers of the field? In my opinion, there 
are three lines of research that will dominate the field in the coming years. 

(a) The local comprehensive simulations will make contact to the global simulations 
of the deep convection zone carried out with anelastic codes. This will require wider and 
deeper boxes and the inclusion of the effects of sphericity and rotation. By combining 
the results (or even the codes), a comprehensive understanding of solar convection and 
its surface manifestations will hopefully be achieved. This is the more important con- 
sidering recent heliose i smic r esults that suggest an inconsistency with global simulations 
( Hanasoge et al.ll2010l [20121 iGizon fe Birchlboil) . 

(&) An second important goal is to extend the simulation domain into the chromo- 
sphere (and corona) with a proper treatment of the radiation processes. In the chromo- 
sphere, the energy balance is dominated by radiative transfer in a few strong spectral 
lines, all of which are formed in NLTE. In addition, ionization and recombination can 
no longer be treated as local equilibrium processes. A consistent treatment of these ef- 
fects within the framework of a 3D MHD simulation provides a c onsiderable challenge , 



which requires a significan t amount of development work (e.g., iGudiksen et al.M2011 



Carlsson fc Leenaartsll2012h . This also i nvolves the treatment of the Hall effect and am- 



bipolar diffusion in non- ideal MHD fe.g- ICheung fc Cameronl2012t iMartmez-Svkora et al 
l2012HKhomenko fc Colladosll2012l ). 

(c) Finally, the comprehensive simulations have reached a degree of maturity and 
validation by solar observations that affords their application to other cool stars. 3D 
hydrodynamic simulations of this kind have b e en carried out already for some time (e.g ., 
Nordlund fc Dravind Il990l : ICollet et all 120071 : iLudwig et alJliooit iFrevtag et all 120121 ). 



Recently, the first magneto-convection sim ulations in the su r face layers of other coo l 



stars than the Sun have become available (jBeeck et alj 120111 : IWedemever et al.l 120121 ) . 



Such simulations will provide a crucial tool for testing and calibrating methods to infer 
stellar magnetic fields from observational spectral and spectro-polarimetric data. 
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